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Abstract 

We compute the exact two-loop matching coefficients for the strong coupHng 
constant ag and the bottom-quark mass mf, within the Minimal Supersymmetric 
Standard Model (MSSM), taking into account 0{al) contributions from Supersym- 
metric Quantum Chromodynamics (SQCD). We find that the explicit mass pattern 
of the supersymmetric particles has a significant impact on the predictions of as and 
rub at high energies. Further on, the three-loop corrections exceed the uncertainty 
due to the current experimental accuracy. In case of the the running bottom-quark 
mass, they can reach in the large tan (] regime up to 30% of the tree-level value. 

PACS numbers: ll.30.Pb, 12.38.-t, 12.38.Bx, 12.10.Kt 



1 Introduction 

Supersymmetry (SUSY) is currently believed to play an important role in physics beyond 
the Standard Model (SM). A compelling argument in favour of SUSY is the particle 
content of the MSSM, that leads in a natural way to the unification of the three gauge 
couplings at a high energy scale /i ~ 10^^ GeV, in agreement with predictions of Grand 
Unification Theories (GUT) [1-3]. 

Apart from the gauge coupling unification, in GUT models based on simple groups such 
as SU{5) [4] or 50(10) [5], also the bottom {yf,) and tau (ur) Yukawa couplings unify at 
the GUT scale. For some models based on 5*0(10) (or larger groups) even the unification 
of the bottom, tau and top (yt) Yukawa couplings is predicted. However, the condition 
of Yukawa coupling unification can be fulfilled within the MSSM only for two regions of 
tan P, the ratio of Higgs field vacuum expectation values: tan /? ^ 1 and tan /? ^ 50 [6-8]. 
A main feature of SUSY models with large tan/3 is that the supersymmetric radiative 
corrections to fermion masses and couplings can be as large as the leading order (LO) 
contributions [7,9,10]. This renders the knowledge of the higher order (HO) corrections 
in perturbation theory mandatory. On the other hand, the unification condition becomes 



very sensitive to the low energy parameters [11]. This property can be exploited to greatly 
constrain the allowed MSSM parameter space. 

With the advent of the CERN Large Hadron Collider (LHC), we will be able to probe 
the realization of SUSY in nature to energy scales of 0{1) TeV. In particular, precision 
measurements and computations will allow to test the low-energy super symmetric rela- 
tions between the particle masses and couplings. It is often argued (for reviews see e.g. 
Refs. [12, 13]) that, from the precise knowledge of the low-energy supersymmetric param- 
eters one can shed light on the origin and mechanism of supersymmetry breaking and 
even on physics at much higher energies, like the GUT scale. The extrapolation of the 
supersymmetric parameters measured at the TeV energy scale to the GUT-scale raises 
inevitably the question of uncertainties involved. Currently, there are four publicly avail- 
able spectrum generating codes [14-17] based on two-loop order MSSM renormalization 
group equations (RGEs) [18-21] subjected to two types of boundary conditions. One set 
of constraints accounts for the weak-scale matching between the MSSM and SM parame- 
ters to one- loop order [22] . The second one allows for the SUSY breaking at the high scale 
according to specific models like minimal supergravity, gauge mediation and anomaly me- 
diation. The approximations within the codes differ by higher order corrections and by 
the treatment of the low-energy threshold corrections. The typical spread of the results is 
within few percents, which does not always meet the experimental accuracy [23]. Along 
the same hne, recent analyses [24, 25] have proven that the three-loop order effects on the 
running of the strong coupling constant and the bottom quark mass may exceed 
those induced by the current experimental accuracy [26,27]. 

Very recently, Ref. [28] computed the two-loop SQCD and top-induced Supersymmetric 
Electroweak (SEW) corrections to the effective bottom Yukawa couplings. The knowledge 
of the two-loop corrections allows predictions for the branching ratios of the MSSM Higgs 
bosons with per-cent level accuracy. 

The aim of this paper is to compute the weak-scale matching relations for the strong 
coupling constant and the bottom quark mass with two-loop accuracy, taking into account 
the exact dependence on the particle masses. This will extend the study of Ref. [24] al- 
lowing phenomenological analyses based on quasi-realistic mass spectra for SUSY par- 
ticles, with three-loop order accuracy. However, we will consider in this paper only 
the 0{a1) corrections from SQCD and postpone the study of the SEW contributions 
0{asy1,0isyliytiy'tyhyb) further investigations. Whereas the SEW corrections to the 
decoupling coefficient of are expected to be negligible, their numerical impact on the 
running bottom quark mass can become as important as those from SQCD. At the one- 
loop order [22], the main contributions to the bottom-quark mass decoupling coefficient 
arise from diagrams containing gluino- and Higgsino-exchange. In the most regions of the 
MSSM parameter space the gluino contribution is the dominant one and can be as large 
as the tree-level bottom quark mass. However, there are domains in the parameter space 
where the corrections due to gluino- and Higgsino-exchange can become of the same order 
and have opposite sign. These regions contain the special points for which the Yukawa 
coupling unification occurs [11,29]. For the MSSM parameters for which the radiative 
corrections to the bottom quark mass are comparable with the LO ones, Ref. [30] proposed 
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a method to resum them to all orders in perturbation theory. A numerical comparison 
with the results of Ref. [30] can be found in Section [51 

The paper is organized as follows: in Section [2] and Section [3] we discuss the theoretical 
framework and the renormalization scheme we use. In Section H] we present the analytical 
one- and two-loop results. The latter ones are displayed in analytical form for three 
simplifying mass hierarchies among the SUSY particles. The numerical effects are studied 
in Section [51 

2 Framework 

As already stated above, the underlying motivation for the running analysis is to re- 
late physical parameters measured at the electroweak scale with the Lagrange parame- 
ters at the GUT scale. The running parameters are most conveniently defined in mass- 
independent renormalization schemes such as MS [31] for the SM parameters and DR [32] 
for the MSSM parameters. It is well known that in such "unphysical" renormalization 
schemes the Appelquist-Carazzone decoupling theorem [33] does not hold in its naive 
form. Quantum corrections to low-energy processes contain logarithmically enhanced 
contributions from heavy particles with masses much greater than the energy-scale of the 
process under consideration. An elegant approach to get rid of this unwanted behaviour 
in the MS or DR scheme is to formulate an effective theory (ET) [34, 35] integrating out 
all heavy particles. The coupling constants defined within the ET must be modified in 
order to account for the effects of the heavy fields. They are related to the parameters of 
the full theory by the so-called matching or decoupling relations. 

For moderate mass splittings between SUSY particles, i. e. there are no large logarithms 
in the theory that have to be resummed by means of RGEs, the decoupling of heavy 
particles might be performed in one step [36]. The energy-scale at which the decoupling 
is performed is not fixed by the theory. It is usually chosen to be /i ~ M, where M is a 
typical heavy particle mass. The MSSM parameters at energies ^ M can be determined 
from the knowledge of the corresponding SM parameters and the associated decoupling 
relations. 

The decoupling coefficients for the strong coupling constant and for the light quark 
masses are known in QCD with four- [37, 38] and three-loop [34] accuracy, respectively. 
Due to the presence of many mass scales, their computation within SQCD and SEW 
is quite involved. At one-loop order, they are known exactly [22,39]. At two-loop 
order, the decoupling coefficient for the strong coupling is known only for specific 
mass hierarchies [24]. Recently, the two-loop SQCD corrections for the decoupling 
coefficient of the bottom-quark mass was computed [39]. The focus of this paper is the 
analytical computation of the decoupling relations for the strong coupling constant and 
the bottom-quark mass within SQCD through two-loops using a different method as 
Ref. [39]. The comparison of the results will be discussed in the next Section. 

We consider SQCD with Uf active quark and Ug = rif active squark flavours and 
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rig = 1 gluinos. Furthermore, we assume that ni = 5 quarks are hght (among which 
the bottom quark) and that the top quark and all squarks and the gluino are heavy. 
Integrating out the heavy fields from the full SQCD Lagrangian, we obtain the Lagrange 
density corresponding to the effective QCD with rii light quarks plus non-renormalizable 
interactions. The latter are suppressed by negative powers of the heavy masses and will 
be neglected here. The effective Lagrangian can be written as follows: 



where ipq, G^, denote the light-quark, the gluon and the ghost fields, respectively, niq 
stands for the light quark masses, ^ is the gauge parameter and Qs = V^vraJ is the strong 
coupling. The index marks bare quantities. C^Q^^ is the usual SQCD Lagrangian from 
which all heavy fields have been discarded. As a result the fields, masses and couplings 
associated with light particles have to be rescaled. They are labeled by a prime in Eq. ([T]) 
and are related with the original parameters through decoupling relations: 

gf = cys , rn^' = c^i , e-i = cue - 1) , 

= ^v^; , G';'^ = vciG','' , cO''-^ = ^/ly'^ . (2) 

Refs. [34] showed that the bare decouphng constants (2, C3 , C3 can be derived from the 
quark, the gluon and the ghost propagators, all evaluated at vanishing external momenta. 
As a result, for calculations performed within the framework of Dimensional Regulariza- 
tion/Reduction (DREG/DRED) only the diagrams containing at least one heavy particle 
inside the loops do not vanish. In Fig. [1] are shown sample Feynman diagrams contributing 
to the decoupling coefficients for the strong couphng (a) and the bottom-quark mass (b). 

For the computation of (g one has to use the well-known Ward identities. A convenient 
choice is the relation: 

^0 ^ ^1 /ox 
^9 7o /To' 

S.3 V ^3 

where Ci denotes the decoupling constant for the ghost-gluon vertex. 
The finite decoupling coefficients are obtained upon the renormalization of the bare pa- 
rameters. They are given by 

/• _ ^^0 ^ _ /.N 

where Z'^ and Z'^ correspond to the renormalization constants in the effective theory. Since 
we are interested in the two-loop results for i = g,m, the corresponding renormalization 
constants for SQCD and QCD have to be implemented with the same accuracy. Analytical 
results for them can be found in Refs. [20,35,40]. 
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(b) 

Figure 1: Sample diagrams contributing to (^3, (^3, ^1 and (m with gluons (g), 
ghosts (c), bottom/top quarks {b/t), bottom/top squarks {b/i) and gluinos {g). 



As mentioned above, the decoupling coefficients can be related with vacuum integrals. 
The latter can be recursively reduced to a master- integral [41] using the method of in- 
tegration by parts [42]. Given the large number of diagrams and occurrence of many 
different mass scales, we computed them with the help of an automated setup. The 
Feynman diagrams were generated with QGRAF [43] and further processed with q2e [44]. 
The reduction of various vacuum integrals to the master integral was performed by a self 
written FORM [45] routine. The reduction of topologies with two different massive and 
one massless lines requires a careful treatment. The related master integral can be easily 
derived from its general expression valid for massive lines, given in Ref. [41]. 



3 Regular izat ion and renormalization scheme 

In our setup, we used the squark mass eigenstates and their mixing angles as input pa- 
rameters. For convenience, we give below the relations between them and the parameters 
of the MSSM Lagrangian. 

The squark mass eigenstates qi^2 ^-nd their mass eigenvalues obtained by diago- 

nalizing the mass matrix 



5 



where we used the notation 



Xg = Aq- jl 



tan /3 , for down-type quarks 
cot [3 , for up-type quarks 



Ml = Ml + ml + MUll-Qqsl,)cos2(3. 



Q 

Ml = Ml^^ + ml + MlQgsl,cos2p. (6) 

Here niq, /| and Qq are the mass, isospin and electric charge of the quark q, respectively, 
and sw = sin6'iy. The parameters Mq and ^ are the soft supersymmetry breaking 
masses, Aq is a trilinear coupling and n is the Higgs-Higgsino bilinear coupling. 
The squark mass eigenvalues are defined through the unitary transformation 



™l, "l-u.A^.ut „,in, sine. 



J -Tl,M,n\, with 7J,= ^J^ ) . (7) 



and the squark mixing angle through 

sin 29, 



2mqXq 



Since we consider the two-loop 0{al) corrections, only the one-loop 0{as) counter- 
terms for the input parameters are required. We have chosen the DR scheme to renormahze 
the strong coupling constant, and the on-shell scheme for the masses of the heavy particles, 
i.e. the gluino, squarks and top quark. The corresponding one-loop renormalization 
constants are known analytically (see, e.g. , Ref. [22]). For the computation of the two- 
loop Cmt, also the one-loop counterterm for the sbottom mixing angle 6b is required. We 
adopted the on-shell renormalization prescription as defined in Ref. [46] 

ReSr (m? ) + ReSr (m? ) 
2[mi -ml) 

where S^^^ is the non-diagonal on-shell sbottom self-energy. 

As we neglect the bottom-quark mass w.r.t. heavy particle masses, an explicit de- 
pendence of the radiative corrections on m^ can occur only through bottom Yukawa 
couplings. In order to avoid the occurrence of large logarithms of the form log(/i^/m^) 
with /i ~ M, we have renormalized the bottom quark mass in the DR scheme. In this way, 
the large logarithms are absorbed into the running mass and the higher order corrections 
are maintained small. 

The renormalization prescription for the trilinear coupling A^ is fixed by the tree-level 
relation Eq. ([8]), while the parameters /i and tan (3 do not acquire 0{as) corrections to 
the one-loop level. Generically, the counterterm for A^ can be expressed as 

5m\ —5m? ^rn^\ "^r —ml 

5Ab= I 2cos29b59b + sm2eb ^ - sm29b — - ^, (10) 

mi — mi mi. j 2mh 

bl 02 / 
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where Sruf, and Sm- are the counterterms corresponding to bottom-quark and squark 
masses, respectively. Due to the use of different renormahzation prescriptions for the 
bottom/sbottom masses and mixing angle, the parameter Ah is renormalized in a mixed 
scheme. 

For the regularization of ultra-violet divergencies, we have implemented DRED with 
the help of the so-called e-scalars [47] . In softly broken SUSY theories, as it is the case of 
MSSM/SQCD, they get a radiatively induced mass. We choose to renormalize their mass 
in the on-shell scheme, requiring that the renormalized mass is equal to zero. 
There are also other approaches available in the literature. We want to mention the one 
proposed in Ref. [48], where the e-scalars are treated as massive particles. This approach 
is known in the literature as the DR' scheme. In this case, the e-scalars have to be 
decoupled together with the heavy particles of the theory [39]. The advantage of this 
method is that it directly relates SQCD parameters regularized within DRED with those 
of QCD regularized within DREG, which are known from experiments. The price of this 
"shortcut" is on one hand, that additional diagrams containing the e-scalars as massive 
particles occur in the calculation of the decoupling coefficients. On the other hand, the 
contributions originating from the change of regularization scheme and those from the 
decoupling of heavy particles are not distinguishable anymore. 

In our approach the change of regularization scheme from DRED to DREG has to be 
performed explicitly within QCD [19,20,49]. For our purposes, the two-loop conversion 
relations for the strong coupling constant and the bottom-quark mass are required. The 
complication that arises at this stage is the occurrence of the evanescent coupling a^. 
associated with the e-scalar-quark-quark vertex. This has to be distinguished from the 
gauge coupling within non-supersymmetric theories (e.g. QCD). However, in SQCD they 
are identical with the gauge couplings, as required by SUSY. Using the ET procedure, we 
can relate them with the strong coupling within the full theory with the help of decoupling 
relations similar with those introduced in Eq. ([2]) 

«e = Cette = Ce«s • (H) 

Following the method described above, one can calculate evaluating the e-scalar and 
quark-propagators and the e-scalar-quark-quark vertex at vanishing external momentum. 
In Ref. [25], its one-loop expression was computed under the simplifying assumption of a 
degenerate SUSY-mass spectrum. In principle, for our numerical analyses, that rely on 
solving a system of coupled differential equations involving also the evanescent coupling 
Oe, even the two- loop order corrections to Ce are needed. However, from the explicit 
calculation it turned out that the numerical effects induced by the two-loop corrections 
to Ce are below the per-mille level. For simplicity, we do not display the corresponding 
two-loop results in the following. The analytical formulae are available upon request from 
the authors. 

The method outlined here and the one introduced in Ref. [39] for the calculation of the 
decoupling coefficient of the bottom-quark mass are equivalent. This has to be understood 
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in the usual sense, that the predictions for physical observables made in one scheme at a 
given order in perturbation theory can be translated to the other scheme through redefini- 
tions of masses and couplings. We have exphcitly checked implementing additionaUy the 
method of Ref. [39] in our setup the equivalence property for the decoupling coefficient 
of the bottom-quark mass through two-loop order. Apart from the obvious rescaling 
of the strong coupling and the bottom quark mass, also the sbottom masses have to be 
modified [48] 

Here is the Casimir invariant in the fundamental representation and denotes the 
mass of the e-scalars. The indices DR and DR', respectively, specify the regularization 
scheme. We also compared numerically the results for the two-loop Cm;, obtained using 
our method with the ones depicted in Fig. 3 of Ref. [39] and found very good agreement. 



4 Analytical results 



4.1 One- loop results 

The exact one-loop results for the decoupling coefficients of the strong coupling constant 
and bottom-quark mass (m can be found in Refs. [22,39]. The analytical formula for 
is new and we give it below up to 0{e). 



(SQCD) 



TT 



- Ic^L, E ^^^^ 
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bi b2 ^ bi 9' 

In Eqs. f lT3l) . f|T4|) . and f|T5|) . we have adopted the abbreviations 



(15) 



Li = ln — , ^ G {t,5f,gi,2,&i,2} , (16) 

where qi{\)i) denote the super-partners of the quark g(6). For our purposes, the special 
case Ce,g=6 is of interest. 

The colour factors are defined in case of a gauge group SU [N) as follows 

- 1 1 
Cf = Ca = N, Tf = -. (17) 

Furthermore, we used the notation ^(2) = tt^/Q and introduced the label "SQCD" to 
specify that the strong coupling has to be evaluated within the full theory, i.e. the SQCD 
with Uf = Ug = 6 active flavours. 

The presence of the terms proportional with the parameter Xb is a manifestation of the 
supersymmetry breaking. They are generated by the Yukawa interaction between left- and 
right-handed bottom squarks and the CP-neutral Higgs flelds. Their contribution to the 
decoupling coefficient of the bottom-quark mass can be related through the Low Energy 
Theorem [50] to the decay rate of the Higgs boson to bb pairs. To one-loop order, the 
Xft-term of Eq. f|T5|) coincides with the SQCD corrections to the decay rate (p ^ bb [51]. 
To higher orders, the relation between the two parameters becomes more involved. 
The Yukawa-coupling induced contributions attracted a lot of attention in the past years, 
due to the fact that they are the dominant corrections for large values of tan (3. They may 
in general become comparable with the tree-level bottom-quark mass. The resummation 
of the one- loop corrections was performed in Ref . [30] . 

4.2 Two- loop calculation 

The analytical two-loop results for the decoupling coefficients are too lengthy 
to be displayed here. They are available in MATHEMATICA format from 
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http://www-ttp.particle.uni-karlsruhe.de/Progdata/ttp08-25 Instead, we present the 
two-loop results for three special cases of the hierarchy among the heavy particle masses. 
Before displaying the analytical results, let us notice the absence of contributions of the 
form a'^X^ to Cmj, in accordance with Refs. [30,51]. They are suppressed by a factor 
rub/M, that we neglect in the ET formalism. 

4.2.1 Scenario A 

We consider first the case of all supersymmetric particles having masses of the same order 
of magnitude and being much heavier than the top-quark 
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M 
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> mt 
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Cs ) Wi, functions of the supersymmetric mass M and the top-quark mass rrit, the soft 



SUSY breaking parameters Xg, q = b,t, the strong coupling constant in the full theory 
^(SQCD) decoupling scale fi. The superscript (5) indicates that the parameters 

are defined in QCD with n; = 5 light quarks. In addition to the notations introduced in 
Eq. (fT6|l . the following abbreviation will be used 
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The two-loop result for the decoupling coefficient of as in case of a degenerate SUSY mass 
spectrum is known [24], however we give it here for completeness 
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Let us point out that, according to Eq. ([8]) the assumption of degenerate top-squark 
masses can be materiahzed only if Xt —>■ 0, due to the heavy top-quark mass. We 
display, however, for completeness the full result. Further on, the hypothesis of equal top- 
and bottom-squark masses is inconsistent with the SU{2) invariance of the t/b isodublet 
imposed in models like mSUGRA. 

4.2.2 Scenario B 

In the following, we discuss the possibility that the gluino is the heaviest supersymmetric 
particle and the squarks have equal masses, much heavier than that of the top-quark 
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The two-loop results read 
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4.2.3 Scenario C 

Finally, we make the assumption that all squark masses are degenerate and are much 
heavier than the gluino and top masses 
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We displayed in the previous expressions only the first three terms of the Taylor 
expansions in the mass ratios. To get an idea about the convergence of the perturbative 
series we fix the following input parameters: mt — 172 A GeV, as^^*^^^ = 0.120, = 
500 GeV, Xg = -4000 GeV, Xt = -400 GeV and let m-g vary. Even for m~g/mg = 0.5, and 
2 the approximations given above agree with the exact results with an accuracy better 
than 1%. For the case of degenerate SUSY masses, i.e m^/m^ = 1 the accuracy is even 
below the per-mille level. 



5 Numerical results 

In this Section we discuss the numerical impact of the two-loop calculations we presented 
above. A first phenomenological application is the prediction of the strong coupling 
and the running bottom-quark mass at high-energy scales like M = 1 TeV or /xqut = 
10^^ GeV, starting from their low-energy values determined experimentally. 

For the energy evolution of the two parameters, we follow the method proposed in 
Ref. [25]: first, we compute ai^\iJ,dec) and m^\ndec) from af\Mz) and mf\mb), re- 
spectively, using the corresponding i-loop SM RGBs [35]. Here fi^cc denotes the energy 
scale at which the heavy particles are supposed to become active, i.e. the scale where 
the matching between the SM and the MSSM is performed. As pointed out in previous 
works [36], one can avoid part of the complications related with the occurrence of the 
evanescent couplings, performing the change of the regularization scheme from the MS 
to the DR scheme at the same scale. Nevertheless, one cannot avoid the occurrence of 
the evanescent coupling ttg in the MS-DR relation for the bottom-quark mass. It has 
to be determined iteratively from the knowledge of the strong coupling at the matching 
scale. For consistency, the i-loop running parameters have to be folded with (i — l)-loop 
conversion and decoupling relations. Above the decoupling scale, the energy dependence 
of the running parameters is governed by the i-loop MSSM RGEs [36]. Wc solved numer- 
ically the system of coupled differential equations arising from the two sets of RGEs, and 
implemented this procedure for i = 1, 2, 3. 

The decouphng scale is not a physical parameter and cannot be predicted by the 
theory. It is usually chosen to be of the order of the heavy particle mass in order to 
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circumvent the appearance of large logarithms. At fixed order perturbation theory, it 
is expected that the relations between the running parameters evaluated at high-energy 
scales and their low-energy values become less sensitive to the choice of /idee once higher 
order radiative corrections are considered. The dependence on the precise value of the 
decoupling scale is interpreted as a measure of the unknown higher order corrections. We 
discuss the scale dependence of as(/iGUT) and mb(yUGUT) in Fig. [2]and Fig. [3l respectively. 

For the SM parameters we used a,(M^) = 0.1189 [26], where Mz = 91.1876 GeV [52], 
mbif^b) = 4.164 GeV [27], with /if, = rubif^b), and Mt = 172.4 GeV [53]. For the SQCD 
parameters, we implemented their values for the SPSla' scenario [13]: rrig = 607.1 GeV, 
TUf^ = 366.5 GeV, m^-^ = 585.5 GeV, m^^ = 506.3 GeV, m^^ = 545.7 GeV, A™(1 TeV) = 
-565.1 GeV, Ap(l TeV) = -943.4 GeV, /i = 396.0 GeV, and tan/5 = 10.0 . 

For the calculation of (mi, to two-loop accuracy, the DR parameter Ab has to be con- 
verted to the renormalization scheme we used here0 For the accuracy level we are con- 
sidering, the one-loop conversion relation is required 

Af'''"^ = Af^ + AAb, where = - "'^^ (25) 

The counterterms SAl were defined in Eq. ( |T0|) and the superscript i indicates the renor- 
malization scheme. The shift AAb is a finite quantity as it can be explicitly checked. It 
depends in turn on the running bottom quark mass in the MSSM. We use an iterative 
method and choose the running bottom quark mass in the SM as the initial parameter. 
A stable solution is obtained after few iterations. In addition, the energy evolution of the 
parameter Ab has to be taken into account. We use here the one- loop RGE, that can be 
derived from Eq. ffTOj) . 

The dependence on the decoupling scale for as(/iGUT) is displayed in Fig. [2l The dot- 
ted, dashed and solid lines denote the one-, two-, and three-loop running, where the corre- 
sponding exact results for the decoupling coefficients have been implemented. One can see 
the improved stability of the three-loop results w.r.t. the decoupling-scale variation. The 
uncertainty induced by the current experimental accuracy on as{Mz), Sas = 0.001 [26], 
is indicated by the hatched band. 

In order to get an idea of the effects induced by the SUSY mass parameters on q;s(/Ugut), 
we show through the dash-dotted line the three-loop results if the SUSY parameters 
corresponding to the Snowmass Point SPS2 [55] are adopted. Their explicit values 
are: m-g = 784.4 GeV, m^-^ = 1003.9 GeV, = 1307.4 GeV, m^^ = 1296.6 GeV, 
= 1520.1 GeV, and tan/5 = 10.0. The curves induced by the other benchmark points 
SPSi, with i = 3, 4, . . . , 9 would lie between the two curves displayed here. One clearly 
notices the great impact of the SUSY-mass pattern on the predicted value of the strong 
coupling at high energies. Accordingly, for precision studies the explicit mass pattern of 
heavy particles must be taken into account. 

In Fig. [3] the scale dependence for mb(/iGUT) is shown. The fine-dotted, dashed and 
solid lines correspond to the exact one-, two-, and three-loop running obtained in the 
SPSla' scenario. As explained above, the energy evolution of the running parameters 

^See Ref. [54] for a detailed discussion. 
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Figure 2: ^^(/iGUT) as a function of /idee- Dotted, dashed and solid lines denote the 
one-, two-, and three-loop contributions, respectively, obtained by using for the input 
parameters their values for the SPSla' benchmark point. The dash-dotted line shows the 
three-loop running corresponding to the SPS2 point. 



have to be combined with the appropriate matching conditions between the low- and 
high-energy regimes. More explicitly, in case of mj, we determine its value within SQCD 
at the energy-scale /^dec through the relation 



f^°(/xdec) = ^^##, Where = . (26) 



Here SCl^'^^ denotes the contributions proportional with tan/5 and ^Cmjf* remaining 
corrections. For simplicity, we do not show in Eq. (!26|) the explicit dependence on the 
MSSM parameters. For the i-loop running analysis, we take into account the {i — l)-loop 
contribution to the Eq. fl26p . As can be seen from the Figure [31 the three-loop results 
stabilize the scale dependence and reduce further the theoretical uncertainty. 
The dotted line displays the two-loop running bottom-mass, where the contributions 
proportional with tan /5 to the one-loop Cmj, are resummed following the method proposed 
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1.3 ^ 




Figure 3: rribificvT) as a function of /idee- The fine-dotted, dashed and sohd hues denote 
the one-, two-, and three-loop contributions to the running bottom mass for the SPSla' 
benchmark point, respectively. The dotted line displays the two-loop running, where the 
tan/3 enhanced contributions are resummed according to Ref. [30]. The dash-dotted line 
represents the three-loop running corresponding to the SPS2 benchmark point. 



in Refs. [13,30]. Within this approach, the matching condition can be written as 

SQCD/ N _ "^r^(A^dec) , 1 _ ^ ~ ^Cmf'^ ^"""^ 

^fe l/^decj - ^i_ioop, : >l-loop - T , r>.tan/3,l-loop " \^ ' ) 



/•i— loop/ \ >i — loop -. I r/-l 



■nib 



The superscript 1 — loop indicates the order in perturbation theory at which the individ- 
ual contributions are evaluated. The authors of Ref. [30] showed that, for a consistent 
analysis not only the tan /3-enhanced contributions have to be resummed, but also the 
next-to-leading logarithms (NLL) a*"*"^ ln*(/i^/m^). In our approach based on z-loop RGEs 
and (i — l)-loop decoupling coefficients the NLL are implicitly resummed. The very good 
agreement between the two computations can be explained by the fact that at one-loop 
order ^Cm** almost an order of magnitude smaller than 5Cm^^- 

The experimental uncertainty generated hj 6 as = 0.001 [26] corresponds to the wider 
hatched band, and the one due to 6 rrib = 25 MeV [27] to the narrow band. Let us notice 
that the three-loop order effects exceed the uncertainty due to current experimental ac- 
curacy 6a s- 
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Finally, the dash-dotted line shows the three-loop running if the SPS2 scenario is im- 
plemented. The differences between the three-loop order results are mainly due to the 
change of masses of the SUSY particles. 









= 1000 GeV 








0.0929 
2.164 


±0.0006|5„4M;,) 

±0.017\sas(Mz) 


±0.015\sm.t{mb) 


-0.003|sPS2 
+0.12 SPS2 


±0.0001|th 
±0.01|th 




f^icn — f^GVT 


Ois (/^rcn) 


0.0405 
1.016 


±0.0001|,„4M;,) 

±0.011\sas{Mz) 




±0.0007|sPS2 
+0.077|sPS2 


±0.0001|th 
±0.005|th 



Table 1: Numerical results for the strong coupling and bottom-quark mass for yUren = 
1000 GeV and /Xren = A'-gut, respectively. The experimental inaccuracy is evaluated 
taking into account the present uncertainty on as{Mz) and mb{mh). The effects of the 
SUSY-mass parameters are evaluated w.r.t. the SPS2 benchmark point. The theoretical 
uncertainties due to unknown higher order corrections are estimated from the variation 
with the decoupling scale for 100 GeV < /idee < 1 TeV. 

For quantitative comparison, we give in Table [T] the numerical values for ^^(/iren) 
and mb(/Xren) for /Xren = 1000 GeV and /iren = AiGUT, evaluated with three-loop accuracy. 
For the decoupling scale we choose /idee = 600 GeV as at this scale the difference 
between the two- and three-loop order corrections reaches a minimum. The different 
sources of uncertainties are exphcitly displayed. The theoretical uncertainties due to 
unknown higher order corrections are estimated from the variation of the three-loop 
results modifying the decoupling scale from 100 GeV to 1 TeV. The effects of the 
SUSY-mass parameters are evaluated as the difference between the three-loop results 
corresponding to the benchmark points SPSla' and SPS2. One can easily see that the im- 
pact of the SUSY-mass pattern is at least five times larger than the experimental accuracy. 

As already pointed out in the previous sections the unification of the Yukawa couplings 
is very sensitive to the MSSM parameters. The dependence on the soft SUSY breaking 
parameters is induced in our approximation only through the decoupling coefficients. 
They comprise an explicit dependence through the Xb parameter (in the case of Cm J and 
an implicit one through the squark masses. 

The analytical formulae for (g and Cm given in Section 14.21 are expressed in terms of 
the physical squark masses. Since they are not known experimentally, for the follow- 
ing numerical analyses, we computed them making the assumption that the soft SUSY 
breaking mass parameters defined in the on-shell scheme obey the following relation 
Mg{t) = Mi){M) = Mfj{M) = Af = ±fi = M, where Mgit) is the on-shell input 
parameter in the stop-mass matrix]^ The corresponding input parameter in the sbottom- 

^See Section |3] for definitions and Refs. [56,57] for a comprehensive discussion. 
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mass matrix acquires a finite shift of 0{as) [56]. Upon diagonalization of the squark-mass 
matrices Eq. ([5]), one obtains the on-shell squark masses mg^^. The parameter enter- 
ing through Xfo the one-loop resuhs have to be converted from the on-shell scheme in the 
renormalization scheme we introduced in Section [21 

In order to estimate the phenomeno logical impact, we discuss the difference between 
m6(/iGUT) evaluated using z-loop running and {i — l)-loop decoupling and the one-loop 
result 

. (i) (i-loop) (1-loop) 

TTlb (1-loop) 

ml 

In Fig. m we fix the soft SUSY breaking mass parameters to M = 1 TeV and Mg = 

800 GeV, and study the dependence of A rnjp for i = 2 (dashed) and i = 3 (solid line) 
as a function of tan (3. One can clearly see the abrupt increase/decrease of the two- and 
three-loop order radiative corrections with the increase of tan/5. Whereas the effects of 
the the one-loop decoupling can be as large as 65% for /i < and —37% for fi > 0, the 
two-loop corrections are moderate, reaching at most 10% and —3%, respectively. 

The numerical effects for large values of tan/3 are of special interest for the study of 

Yukawa-coupling unification. In Fig. [5] we show A m^^^ as a function of the soft SUSY 
breaking mass scale for tan /3 = 50. One can see the increase in size of radiative corrections 
for lighter SUSY masses. Again, for n > the bulk of the corrections are comprised in 
the two- loop running mass, while the three- loop order effects sum up to few percent. For 
/i < and light SUSY masses the three-loop contributions can increase the bottom-quark 
mass with up to 30%. 



6 Conclusions 

The knowledge of fundamental parameters at high energies, such as M or /xgut, are 
essential for the reconstruction of the theory beyond the SM. In this paper we presented 
the exact two-loop decoupling coefficients of the strong coupling and the bottom-quark 
mass within the SQCD. Together with the known three-loop order RGBs they allow 
predictions of the two parameters at high energies with three-loop accuracy. This level of 
precision on the theory side is necessary in order to match with the current experimental 
accuracy. The values of the gauge and Yukawa couplings at the unification scale /icuT 
are essential ingredients for the determination of the GUT threshold corrections, which 
in turn are used to identify the underlying GUT model. 

In addition, the dependence on the energy scale at which the supersymmetric particles 
are integrated out, which reflects the size of the unknown higher order corrections, is 
significantly reduced in the case of the three-loop order predictions. 

Furthermore, the approach outlined here accounts for the effects induced by the in- 
dividual mass parameters. They are phenomenologically significant for both parameters 
and exceed the experimental uncertainty by more than a factor five. 
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tan P 

(b) 

Figure 4: Arrib/mb as a function of tan/3 for /i < (a) and /i > (b). The 
soft SUSY breaking mass parameters are fixed to 1 TeV and the gluino mass 
to Mg = 800 GeV. 



21 





22 



The radiative corrections to the running bottom-quark mass are particularly important 
for SUSY models with large values of tan/?. It turns out that for negative values of /i 
the three-loop order contributions can reach up to 30% from the tree-level bottom-quark 
mass. Furthermore, they clearly stabilize the perturbative behaviour. These features 
render them indispensable for studies concerning the Yukawa- coupling unification. 
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